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SUMMARY 

We  analyze  the  PP  phase  at  upper  mantle  distances  (25°  to  60°)  to  quantify  its 
propagation  characteristics  and  potential  for  determining  cmst  and  upper  mantle  velocity 
structure.  Upper  mantle  distance  PP  waveforms  are  very  complex,  involving  interference 
between  a  variety  of  arrivals  traversing  different  depth  ranges  in  the  Earth.  These  arrivals 
include  PP  triplications  produced  by  the  two  major  (and  any  other)  upper  mantle  seismic 
discontinuities;  energy  that  turns  in  the  upper  mantle  and  reflects  from  the  underside  of  the 
Moho  (e.g.  PmP);  and  multiply  reflected  or  converted  phases  in  the  crust  at  the  source, 
receive  and  PP  mid-path  bounce  point  (PP  coupled  PL).  Even  for  laterally  homogeneous 
structures,  con:q>iete  synthetic  seismograms  such  as  provided  by  reflectivity  calculations  are 
required  to  adequately  model  the  full  suite  of  arrivals  that  contribute  to  upper  mantle  PP 
phases.  TTie  complexity  of  PP  waveforms  can  be  exploited  to  extract  crust  and  upper 
mantle  structural  information  for  relatively  homogeneous  paths,  but  the  numerous 
interfering  phases  make  it  difficult  to  attribute  PP  characteristics  to  panicular  structural 
features  without  some  a  priori  information.  PP-P  differential  travel  times  and  PP/P 
amplitude  ratios  provide  constraints  on  average  upper  mantle  velocities  above  tlie  transition 
zone  if  the  general  crustal  properties  along  the  path  are  known.  Often  a  priori  information 
on  crustal  thickness  is  available  and  average  lid  velocity  structure  can  thus  be  determined. 
Large  amplitude  early  PP  coda,  comprised  of  crustal  reflections  and  conversions,  can 
provide  information  on  crustal  thickness  and  must  be  accounted  for  when  identifying 
secondary  upper  mantle  PP  triplication  branches. 

It  is  also  possible  to  place  some  constraints  on  the  velocity  gradients  in  the 
uppermost  mantle.  In  the  presence  of  a  positive  velocity  gradient  below  the  crust, 
whispering  gallery  multiple  underside  reflections  from  the  Moho  can  be  ob.served  at  uppei' 
mantle  distances,  "nie  amplitudes  of  these  Moho  whispering  gallery  arrivals  and  associate! 
underside  reflections  (PmP>  PmPmP>  etc.)  are  sensitive  to  the  lid  velocity  gradient, 
particularly  for  short-period  and  broadband  data  at  distances  from  25°  to  38°.  We 
demonstrate  the  potential  of  accurate  PP  waveform  modeling  for  detemrining  upper  mantle 
structure  by  modeling  long-period  signals  traversing  relatively  'pure'  paths  across  North 
America.  A  continuum  of  P  wave  velocity  models  with  varying  lid  structure  is  appropi  ititc 
for  the  upper  mantle  beneath  North  America.  The  highest  velocities  are  found  beneatli  the 
northern  Canadian  shield,  with  the  lid  becoming  progressively  lower  velocity  beneath  the 
southern  shield,  continental  platform  and  tectonically  active  regions  respectively.  Most  of 
the  heterogeneity  in  P  velocity  structure  is  concentrated  above  a  depth  of  250  km,  alihoiigli 
some  data  are  consistent  with  models  that  are  relatively  slow  at  greater  depths.  Refinement 
of  these  models,  and  full  exploitation  of  the  PP  phase  will  require  extensive  brotidband  d:it:i 


sets  and  two-  and  three-dimensional  synthetic  capabilities  that  can  handle  the  intrinsic 
complexity  of  PP  in  the  presence  oflatcr.il  heterogeneity. 


I  INTRODlJCnON 


Calculation  of  synthetic  seismograms  has  contributed  significantly  to  our 
knowledge  of  the  Eanh's  velocity  structure  by  allowing  waveforms  from  many  diverse 
phases  to  be  successfully  modeled  or  inverted  for  deep  structure.  A  variety  of  techniciiies 
are  now  available  for  the  computation  of  theoretical  body  wave  seismograms  in  realistic 
Earth  models  and  the  choice  of  method  for  a  particular  application  is  usually  determined 
through  a  compromise  between  the  desired  accuracy  and  the  time  and  expense  required  to 
perform  the  calculation.  The  WKBJ  algorithm  in  its  original  form  (Chapman  1978)  is 
exceptionally  fast  but  involves  significant  approximations,  reducing  its  suitability  for  some 
applications.  'liie  Cagniard-de  Hoop  method  (e.g.,  ilelmberger  1973,  1983)  has  fewer 
limitations  and  provides  superb  insight  into  the  seismogram  construction,  but  like  the 
WKBJ  method  it  requires  specification  of  the  ray.set  to  be  considered.  3'he  reflectivity 
technique  (Fuchs  &  Muller  1971)  employs  propagator  matrices  to  account  for  all 
conversions  and  reflections  in  a  stratified  medium,  and  is  thus  quite  accurate  and  complete, 
but  the  method  is  relatively  costly  and  it  is  often  difficult  to  gain  insight  into  the  nature  of 
individual  arrivals.  These  procedures  have  been  extensively  compared  and  their  respective 
limitations  are  well  understood  in  general  (e.g.,  Burdick  &  Orcutt  1979;  Chapman  and 
Orcutt  1985).  Many  additional  procedures  and  hybrid  methods  have  been  developed  for 
application  in  both  layered  and  laterally  heterogeneous  structures.  Richards  (1985) 
provides  a  review  of  many  applications  of  synthetic  seismogram  modeling  for  determining 
structure  of  the  crust,  mantle,  and  core. 

The  principal  phases  used  for  characterizing  upper  mantle  structure  have  been  direct 
P  and  S  pha.ses  that  interact  with  the  transition  zone  velocity  discontinuities  to  produce 
upper  mantle  triplications,  reflections  and  conversions.  In  order  to  expand  the  spatial 
coverage  of  upper  mantle  structure,  seisniolcjgists  have  begun  to  model  later  phti.ses, 
involving  multiple  surface  and/or  core  reflections  along  their  path.  The  application  of 
WKBJ  and  Cagniard-de  Hoope  synthetic  modeling  to  multiple  bounce  transverse  shear 
waves  (SS  and  SSS)  was  instrumental  iq  establishing  the  magnitude  and  depth  extent  of 
lateral  variations  in  the  upper  mantle  shear  wave  velocity  structure  beneath  North  America 
(Grand  &  Helmberger  1984a)  and  the  North  Atlantic  Ocean  (Grand  &  Helmberger  1984b). 
The  success  of  this  SS  waveform  modeling  in  revealing  the  .shear  wave  structure  of  the 
upper  mantle  beneath  North  America  led  LeFevre  &  Helmberger  (1989)  to  include  PP 
phases  at  upper  mantle  distances  with  direct  P  and  crustal  P^/  pha.ses  in  developing  an 
upper  mantle  P  wave  model  (S25)  appropriate  for  the  Canadian  Shield.  WKBJ  and 
Cagniard-de  Hoope  synthetics  for  P  and  Pnl  waves  for  model  S25  fit  the  observations  very 
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satisfactorily.  However,  simplified  WKBJ  synthetics,  involving  only  'piimary'  ray  sets 
like  those  used  to  model  the  direct  P  phases,  do  not  match  many  of  the  PP  data  nearly  as 
well.  This  appears  to  result  from  incomplete  synthesis  of  the  PP  arrival  rather  than 
inadequacy  of  model  S25.  The  contrasting  success  of  simple  WKBJ  modeling  of 
transverse  SS  wavefoitns  (Grand  &  Helmberger,  1984a)  further  indicates  that  the  P-SV 
system  for  multiple  surface  adlections  is  more  complicated  to  mixlel. 

In  this  paper  we  perform  a  detailed  analysis  of  the  PP  phase  at  upper  maiiiic 
distances  (25"  to  60°;  twice  the  corresponding  upper  mantle  distance  range  for  diieci  P)  to 
better  understand  its  propagation  and  the  requirements  for  accurate  synthetic  modeling  of 
the  phase.  We  then  use  this  insight  to  evaluate  the  usefulness  of  PP  modeling  in 
determining  upper  mantle ‘structure.  We  find  that  the  previous  difficulty  in  matching 
observed  PP  waveforms  is  largely  due  to  inadequacies  in  the  simplified  WKBJ 
computations  that  were  made.  Complete  reflectivity  synthetic  seismograms,  or  more 
extensive  WKBJ  synthetics,  can  ade<|uatcly  acc«>nnt  for  the  complexiiy  ol  tipper  mantle  IM' 
phases,  and  provide  a  significantly  iinproved  fit  to  the  data  that  enables  fniiher  recovery  of 
earth  structure  information.  We  present  a  variety  of  synthetic  examples  to  elucidate  the 
characteristics  of  the  PP  phase  and  establish  its  potential  utility  for  interrogation  of  upper 
mantle  velocity  structure.  We  conclude  with  preliminary  modeling  of  PP  phases  traversing 
the  upper  mantle  beneath  North  America  in  which  we  resolve  lateral  variations  in  upper 
iruintlc  P  wave  velocities  beneath  the  continent. 

2  SYNTHETIC  PP  WAVEFORMS  AT  UPPER  MANTLE  DISTANCES 

Grand  &  Helmberger  (1984a)  and  Lefevre  &  Helmberger  (1989)  have 
demonstrated  that  multiple  surface  reflections  (PP,  SS,  SSS  etc.)  that  turn  in  the  upper 
mantle  transition  zone  can  be  used  to  greatly  extend  our  coverage  of  upper  mantle  structure. 
The  jirimary  advantages  of  these  phases  aic  that  the  tip|KT  mantle  tiiphcations  of  I’l’.iiul  S.s 
arc  at  twice  tlie  disttuice  (25°  to  60°)  and  have  twice  the  lime  separation  between  triplication 
branches  compared  to  the  associated  triplications  of  the  direct  P  and  S  phases  (obsers  ed  at 
distances  of  12.5°  to  30°).  As  a  result,  the.se  phases  have  been  valuable  in  determining 
velocity  models  for  regions  like  the  Canadian  .shield  and  East  Pacific  Rise,  where  only 
sparse  direct  pha.se  coverag.e  :it  upper  mantle  tii.siances  is  available.  The  primary  focus  of 
previous  re.search  has  Ikvii  on  identifying  the  upper  mantle  triplication  (and  ijuintuplicalion) 
arrivals  in  the  multiple  surface  reflections.  1-igure  I  shows  typical  veriictil  component 
long-period  WWSSN  and  Canadian  Seismic  Network  recordings  of  P  and  PP  waveforms 
from  the  October  18,  1967  Arctic  event.  These  data  were  modeled  by  LeFevre  & 
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Helmberger  (1989)  in  the  derivation  of  their  Canadian  shield  P  wave  model,  S25.  These 
waveforms  are  at  distances  where  direct  P  waves  have  bottoming  points  in  the  lower 
mantle,  and  any  information  about  the  upper  mantle  structure  along  the  path  must  be 
derived  from  the  shallower  turning  energy  associated  with  the  PP  phase. 

LeFevre  &  Helmberger  (1989)  calculated  synthetic  seismograms  for  the  1967  Arctic 
event  using  the  WKBJ  technique  (Chapman  1978).  In  Figure  1,  we  show  very  similar 
synthetics  for  model  S25  calculated  in  the  same  manner.  Inspection  of  the  observed 
waveforms  and  the  WKBJ  synthetics  reveals  that  although  the  PP-P  differential  travel 
times  are  generally  matched  by  this  model,  the  synthetic  PP  waveforms  do  not  clo.seiy 
resemble  the  data.  The  triplication  arrivals  predicted  by  the  WKBJ  synthetics  for  model 
S25  arrive  in  a  narrow  time  window  at  distances  from  40°  to  45°  (see  stations  SFA  and 
OTT  in  Figure  1 )  resulting  Tn  large  amplitude,  relatively  simple  PP  waveforms.  Although 
the  data  do  exhibit  fairly  simple  waveforms  at  these  stations,  the  polarities  of  the  largest 
amplitude  arrivals  are  actually  reversed  relative  to  the  synthetics,  and  one  would  be  hard- 
pressed  to  reliably  identify  the  triplication  features.  At  larger  distances,  the  EF  branch 
(energy  turning  below  670  km)  and  the  AB  branch  (energy  turning  above  4(X)  km)  of  the 
upper  mantle  triplications  separate.  This  separation  is  apparent  in  the  data,  however,  there 
are  later  arrivals  with  comparable  amplitude  that  are  not  predicted  by  the  WKBJ  synthetics 
(e.g.,  OGD  and  GEO).  This  degree  of  waveform  mismatch  between  data  and  synthetics 
makes  it  difficult  to  identify  the  onset  of  the  emergent  PP  p’.iase  and  its  triplicated  arrivals. 
The  quality  of  the  WKBJ  synthetic  fit  to  the  PP  data  alone  would  provide  little  confidence 
in  the  details  of  model  S25;  LeFevre  &  Helmberger  (1989)  relied  primarily  on  direct  P  and 
P„l  observations  in  the  derivail.  ’  of  their  model,  with  the  PP-P  timing  being  u.sed 
essentially  (and  successfully)  as  a  consistency  check. 

Synthetic  waveforms  computed  with  the  reflectivity  method  for  model  S25  differ 
dramatically  from  the  simple  WKBJ  synthetics,  as  shown  in  Figure  1,  and  provide  much 
better  fits  to  the  PP  data.  Although  the  waveform  modeling  requires  further  improvement, 
the  polarity  of  the  PP  waveforms  at  SFA  and  OTT  is  now  correctly  predicted,  and  the  large 
amplitude  arrivals  following  the  AB  branch  at  OGD,  GEO,  FLO  and  OXF  are -also 
accounted  for.  The  comparison  at  OTT  is  particularly  striking.  The  same  triplications  are 
present  in  the  two  synthetics,  but  the  reflectivity  synthetic  predicts  so  much  additional 
energy  arriving  in  the  overall  PP  waveform  that  even  identifying  the  triplication  branches  is 
very  difficult.  However,  the  fact  that  the  reflectivity  synthetic  re.sembles  the  data  quite 
closely  holds  out  hope  for  using  the  PP  waveforms  to  study  upper  mantle  structure.  The 
reflectivity  synthetics  also  produce  a  coda  for  direct  P  which  was  not  included  in  the 
WKBJ  synthetics,  and  which  does  not  match  the  data  veiy-  well.  While  the  reflectivity 
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procedure  is  restricted  to  laterally  homogeneous  layers,  and  the  crustal  structure  may  in  fact 
vaiy  along  the  path,  it  is  not  surprising  that  model  S25  produces  such  strong  P  coda.  Later 
small  arrivals  between  P  and  PP  in  the  reflectivity  synthetics  are  produced  by  underside 
reflections  from  upper  mantle  discontinuities  that  will  be  discussed  below. 

Grand  &  Helmberger  (19ti4a)  compared  WKBJ  synthetics  with  more  accurate 
Cagniard-de  Hoope  generalized  ray  calculations  (Wiggins  &  Helmberger  1974)  for 
transverse  SS  phases  at  upper  mantle  distances  and  did  not  find  a  comparable  extent  of 
waveform  incompatibility.  Of  course,  both  methodologies  are  ray-based,  and  similar 
truncated  ray  expansions  were  used  in  each  method.  The  observation  that  the  synthetics 
reproduce  much  of  the  observed  waveform  character  does  suggest  that  truncated  ray 
expansions  are  sufficient  for  the  SS  phase.  The  PP  phase  is  more  complicated  than 
transverse  SS  at  upper  mantle  distances  due  to  multiple  P-SV  conversions  from  interfaces 
near  the  source,  receiver,  and  PP  mid-path  bounce  point.  The  WKBJ  synthetics  computed 
by  LeFevre  &  Helmberger  (1989)  and  our  own  shown  in  Figure  1  include  only  the  primarj' 
P  and  triplicated  PP  waves  turning  in  the  mantle  and  their  corresponding  depth  phases. 
Phases  reflected  and/or  converted  from  the  Moho  or  strong  transition  zones  within  the  crust 
are  not  included  in  these  WKBJ  synthetics,  thus  the  synthetics  are  a  partial  calculation  of 
the  complete  layered  Earth  response.  It  is  well-established  that  at  distances  beyond  70°, 
where  the  PP  phase  turns  in  the  lower  mantle,  simplified  primary  ray  expansions  in 
conjunction  with  a  Hilbert  transform  produce  PP  synthetics  that  match  observed  data  ven' 
well  (e.g.,  Lynnes  &  Ruff  1985),  but  this  apparently  is  not  tme  at  upper  mantle  distances. 
We  found  the  degree  of  waveform  error  in  the  truncated  synthetic  calculations  to  be 
surprising,  and  were  thereby  motivated  to  conduct  this  analysis. 

Figure  2  isolates  the  primary'  structural  features  affecting  the  discrepancy  between 
the  primary  ray  and  complete  synthetics  by  comparing  P  and  PP  synthetics  calculated  with 
WKBJ  and  reflectivity  for  simplified  earth  models.  We  make  this  comparison  using  the 
reflectivity  method  rather  than  generalized  rays  so  Uiat  we  do  not  need  to  specify  a  ray.sei 
and  are  ensured  of  including  all  imponant  reverberations.  Tiiese  synthetics,  and  all  others 
discussed  in  this  section  unless  otherwise  stated,  are  vertical  component  seismograms 
computed  for  a  vertical  dip-slip  source  (strike=0o,  dip=90o,  rake=90o)  at  a  depth  of  1 8  km 
and  convolved  with  a  long-period  WWSSN  insmument  response.  The  attenuation  model 
from  PREM  (Dziewonski  &  Anderson  1981)  is  used  in  the  computation  of  the  reflectivity 
synthetics,  while  attenuation  is  included  in  the  WKBJ  synthetics  by  convolution  with  an 
anelastic  attenuation  operator  (Futierman  1962)  having  a  t*  value  of  1.0  s.  The  synthetics 
in  Figure  2  were  calculated  using  the  four  different  velocity  models  indicated:  a  smooth 
positive  velocity  gradient  with  no  seismic  discontinuities  and  no  crustal  layer  (grad  no 
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crust),  the  same  model  with  a  low  velocity  crustal  layer  (grad),  a  model  with  increasing 
velocities  with  depth  with  two  upper  mantle  seismic  discontinuities  and  no  crust  (disc  no 
crust),  and  the  same  model  with  a  crust  (disc).  The  correspondence  between  synthetics 
computed  for  models  that  do  not  contain  a  low  velocity  crustal  layer  (grad  no  crust  and  disc 
no  crust)  is  excellent,  demonstrating  that  triplication  effects  are  well  accounted  for  by  the 
WKBJ  method  as  expected.  However,  the  synthetics  computed  for  models  with  a  crustal 
layer  are  quite  different.  The  largest  discrepancy  between  the  synthetics  shown  in  Figure  2 
occurs  in  the  later  part  of  the  P  and  PP  waveforms,  indicating  that  the  extra  energy  in  the 
complete  synthetics  arises  largely  from  multiple  bounces  in  the  crust.  The  waveform 
di.screpancies  are  more  pronounced  for  the  model  with  upper  mantle  discontinuities  becaii.se 
the  multiple  reverberations  also  triplicate,  producing  strong  interference  patterns.  At 
distances  larger  than  45°  this  crustal  complexity  is  manifested  sufficiently  late  in  the 
waveforms  that  the  early  part  of  the  PP  phase  is  not  distorted,  but  this  is  clearly  not  the 
case  at  42°,  where  the  interference  almost  reverses  the  first  strong  PP  upswing.  This  can 
easily  lead  to  confusion  when  picking  arrivals  in  the  data. 

In  addition  to  crustal  reverberations  that  dominate  the  later  part  of  the  PP 
waveforms,  the  Moho  will  also  produce  PP  arrivals  that  reflect  from  its  underside  rather 
than  from  the  free  surface  (PmP)  and  arrive  prior  to  PP.  Figure  3  shows  the  influence  of 
these  arrivals  on  PP  waveforms  as  a  function  of  propagation  distance.  Synthetic 
seismograms  are  calculated  with  the  negative  gradient  velocity  model  illustrated  in  Figure 
4.  Waveforms  in  the  top  three  rows  are  WKBJ  synthetics  that  include  the  following  rays: 
direct  P,  pP,  sP  and  PP  precursors  reflecting  up  to  three  times  from  the  bottom  side  of  the 
Moho  (PmP,  Pm  PmP,  PmPmPmP)i  direct  P,  pP,  sP,  PP,  pPP  and  sPP  (PP);  and  a 
combination  of  the  previous  2  raysets  (PmP+PP).  The  Moho  reflected  PP  energy  precedes 
the  PP  phase  out  to  distances  of  about  40°,  producing  waveforms  (PmP+PP)  that  begin 
earlier  and  differ  significantly  from  synthetics  computed  with  the  PP  phase  alone  (PP). 
Beyond  40®,  the  Moho  reflected  phases  have  reduced  amplitude  and  arrive  simultaneously 
with  PP,  changing  it.s  waveform  less  dramatically.  The  bottom  row  in  Figure  3  shows 
reflectivity  synthetics  which  inherently  include  all  of  the  reverberations.  Comparison 
between  WKBJ  and  reflectivity  synthetics  in  the  distance  range  where  Moho  reflected 
phases  contribute  to  the  PP  waveforms  indicates  that  these  arrivals  are  indeed  important, 
and  are  partially  responsible  for  the  incompatibility  between  the  synthetic  seismograms. 
The  arrival  between  P  and  PP  visible  in  all  WKBJ  synthetics  but  absent  from  the 
reflectivity  synthetics  is  the  B  branch  of  the  P  wave  triplication.  The  overestimation  of 
WKBJ  synthetic  amplitudes  in  the  shadow  of  triplications  results  from  inadequacies  of 
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Fresnel  diffraction  theory  at  the  slowness  of  grazing  rays  and  is  well-established  (e.g.. 
Chapman  &  Orcutt  1985). 

2.1  Whispering  Gallery  Modes 

For  a  positive  velocity  gradient  in  the  uppermost  mantle  beneath  the  crust,  Menke  & 
Richards  (1980)  showed  that  the  Moho  can  act  as  an  efficient  'whispering  galler/ 
interface,  causing  energy  to  propagate  along  it.  We  follow  the  convention  of  identifying 
phases  as  'whispering  gallery'  arrivals,  PmPfWG),  only  if  they  'hug'  the  boundar>'  at 
which  they  have  underside  reflections.  Phases  turning  deeper  in  the  upper  mantle  and 
reflecting  from  the  underside  of  the  Moho  are  termed  PmP  as  discussed  in  the  previous 
section.  Figure  5  schematically  illustrates  these  phases.  The  generation  of  significant 
PinP(WG)  requites  that  the  impedance  contrast  across  the  Moho  be  large  enough  to  reflect 
the  energy  turned  in  the  positive  gradient.  Menke  &  Richards  (1980)  demonstrated  that  for 
Eanh  models  satisfying  this  criteria,  the  family  of  Moho  whispering  gallery  phases  can 
have  significant  amplitude  relative  to  the  direct  P  wave.  At  upper  mantle  PP  distances 
(250-600),  Moho  whispering  gallery  phases  can  arrive  between  15  and  150  s  after  the 
direct  P  wave,  ahead  of  and  overlapping  the  PP  phase,  and  thus  PmPfWG)  may  also 
contribute  to  the  mismatch  apparent  in  Figure  1  between  WKBJ  and  reflectivity  synthetics. 

Menke  &  Richards  (1980)  found  that  the  velocity  gradient  within  the  lid  is  the  most 
important  parameter  affecting  the  properties  of  PniP(WG).  Therefore  at  certain  distance.^, 
where  the  early  part  of  the  PP  phase  consists  of  whispering  gallery  arrivals,  thorough 
analysis  of  PP  waveforms  may  provide  constraints  on  lid  velocity  gradient.  We  investigate 
this  by  comparing  reflectivity  synthetics  computed  for  velocity  models  having  posiii\e 
gradient,  negative  gradient  and  constant  lid  velocities  (Figure  4).  The  gradients  are 
approximated  by  many  thin  homogeneous  layers  in  the  reflectivity  calculations. 

Rerlectivity  synthetics  for  the  models  in  Figure  4  are  shown  in  Figures  6  (long- 
period)  and  7  (broadband).  For  long-period  waveforms  out  to  a  distance  of  34°  strong, 
impulsive  arrivals  are  apparent  between  P  and  PP  for  the  positive  lid  gradient  synthetic.s, 
with  weaker  energy  arriving  at  the  same  time  in  the  constant  lid  velocity  synthetics,  and 
little  energy  being  present  in  the  waveforms  for  the  negative  lid  gradient  model  (Figure  6). 
The  impulsive  arrivals  in  this  di.stance  range  for  the  positive  lid  gradient  involve  whispering 
gallery  phases  that  turn  in  the  upper  165  km  of  the  mantle.  The  whispering  gallery'  phases 
are  not  observed  for  the  negative  gradient  velocity  m.odel  since  no  energy  turns  in  the  top 
165  km.  Beyond  34®,  energy  bottoming  in  the  deeper  positive  velocity  gradient  below  165 
km  (figure  4)  and  reflected  from  the  boitomside  of  the  Moho  beco.mes  the  dominant  arrival 
precaling  PP  for  all  of  the  velocity  models.  This  phase  corresponds  to  P^P  in  Figure  3. 
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and  the  long-period  signals  beyond  34°  are  not  diagnostic  of  the  differences  in  lid  velocity 
gradient.  At  distances  closer  than  25®,  the  whispering  gallery  arrivals  produced  by  the 
positive  lid  velocity  gradient  model  a'-e  obscured  by  the  relatively  large  secondary  branches 
of  the  direct  P  wave  upper  mantle  triplication  arrivals  and  P  coda.  Therefore  long-period 
waveforms  in  the  distance  range  between  26®  to  34®  are  most  sensitive  to  the  velocity 
gradient  in  the  lid,  but  it  should  be  noted  that  our  models  are  rather  extreme.  The  more 
subtle  long-period  waveform  differences  for  smaller  variations  in  lid  gradient  could  easily 
be  missed  in  real  data  where  lower  signal  to  noise  ratios  are  expected. 

Menke  and  Richards  (1980)  showed  that  whispering  gallery  phases  tend  to  lose 
much  of  their  low-frequency  energy  through  tunneling  in  tlic  lower  velocity  material 
underlying  the  lid,  making  PniP(WG)  more  prominent  at  shorter  periods.  This  effect  is 
apparent  in  broadband  synthetics  for  the  same  models  described  above  (Figure  7).  The 
broadband  synthetics  for  the  positive  lid  velocity  gradient  model  have  PmP(WG)  arrivals 
with  about  ten  times  larger  amplitude  than  the  energy  arriving  at  the  same  time  for  the 
negative  lid  velocity  gradient  model  and  twice  the  amplituae  of  the  energy  in  the  constant 
lid  velocity  synthetics  for  distances  out  to  34°.  At  distances  between  35°  and  38°  the 
negative  lid  velocity  gradient  synthetics  have  larger  amplitudes  in  the  early  part  of  the  PP 
waveforms  than  synthetics  for  the  other  velocity  models  (Figure  7).  This  energy  turns  in 
the  positive  velocity  gradient  above  the  400  km  discontinuity  and  reflects  from  the 
bottomside  of  the  Moho  (Pm?)-  The  energy  that  was  not  turned  in  the  lid  for  the  negative 
lid  gradient  model  strengthens  the  arrivals  in  this  distance  range,  providing  another 
possible  diagnostic  of  the  lid  structure  (note  that  in  these  calculations  all  of  the  models  have 
similar  velocity  gradients  from  200  to  400  km  depth;  see  Figure  4).  Beyond  38®,  little 
PmP  precursory  energy  persists  and  synthetics  for  all  of  the  velocity  models  resemble  one 
another. 

These  calculations  suggest  that  in  the  distance  range  26®  to  38®  the  variation  in 
broadband  PP  waveforms  predicted  for  upper  mantle  models  with  .significantly  different  lid 
gradients  is  certainly  large  enough  to  detect  in  "real  data,  if  homogeneous  paths  are 
traversed.  The  variation  in  lid  gradients  represented  in  Figure  4  is  almost  certainly  larger 
than  we  expect  in  the  Earth;  however,  this  exercise  illustrates  that  analysis  of  the  complete 
PP  phase  can  potentially  differentiate  between  positive  and  negative  upper  mantle  lid 
gradients.  Determination  of  the  lid  gradient  is  of  profound  importance  for  the  high 
frequency  propagation  of  regional  phases  like  P^  (e.g.,  Sereno  &  Given  1990),  which  are 
used  for  earthquake  and  nuclear  explosion  discrimination.  At  larger  distances,  from  38°  to 
60°,  energy  that  arrives  between  P  and  PP  is  mainly  caused  by  underside  reflections  from 
deeper  discontinuities  (PdP).  including  the  400  and  670  km  discontinuities.  Neele  & 
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Snieder  (1991)  demonstrate  that  this  late  P  coda  is  quite  coherent  across  small  arrays,  and 
waveform  modeling  of  the  PdP  phases  can  be  used  to  constrain  the  discontinuities  at  the 
base  of  the  lid  and  in  the  transition  zone. 

22  Crustal  Reverberations 

Although  explicit  inclusion  of  whispering  gallery  modes  and  other  bottomside 
Moho  reflections  in  the  WKBJ  synthetics  improves  comparison  with  reflectivity  synthetics 
for  the  early  part  of  the  PP  waveforms  (Figure  3),  the  high  amplitude  later  part  of  the 
reflectivity  PP  waveforms  is  not  yet  accounted  for.  To  determine  what  phases  contribute  to 
the  large  PP  coda  in  the  reflectivity  calculations  (and  actual  data),  we  added  additional  rays 
to  the  WKBJ  synthetics  (Figure  8).  The  synthetics  in  Figure  8  were  computed  using  a 
velocity  model  with  a  slight  positive  gradient  in  the  upper  2(X)  km  of  the  mantle  overlain  by 
a  6.4km/s,  40  km  thick  crust.  The  synthetics  in  the  top  row  of  Figure  8  (WKBJl)  were 
computed  with  the  WKBJ  procedure  including  only  primary  direct  P  and  PP  phases  and 
their  associated  depth  phases,  pP,  sP,  pPP,  and  sPP.  In  addition  to  this  rayset,  the 
WKBJ2  synthetics  include  rays  turning  in  the  upper  400  km  of  the  mantle  and  bouncing 
once  from  the  underside  of  the  Moho  (this  includes  the  first  whispering  gallery  mode 
which  is  small  due  to  the  weak  velocity  gradient  in  the  lid)  and  all  rays  having  an  extra 
reflection  or  conversion  from  both  the  free-surface  and  the  Moho.  The  bottomside  Moho 
reflected  energy  contributes  to  the  beginning  of  the  PP  waveform  as  discussed  in  the 
previous  section. 

Inclusion  of  first  crustal  multiples  begins  to  build  the  large  amplitudes  observed  in 
the  later  part  of  the  reflectivity  PP  waveforms,  bringing  the  WKBJ  and  reflectivity 
synthetics  into  better  agreement  (Figure  8).  This  is  especially  apparent  at  distances  larger 
than  40°  where  the  waveforms  are  relatively  simple.  The  P  waveforms  remain  unchanged 
since  the  P  wave  crustal  multiples  were  not  included  in  the  WKBJ  calculations  so  that  we 
could  focus  on  the  later  arrivals.  The  agreement  between  WKBJ  and  reflectivity  synthetics 
clearly  improves  with  the  addition  of  rays  that  have  only  one  extra  bounce  in  the  crust. 
Additional  improvement  is  expected  as  more  crustal  reflections  and  conversions  are  added 
until  all  possible  interactions  in  a  layered  Earth  have  been  accounted  for.  But, 
unfortunately,  the  convergence  is  rather  slow  at  distances  less  than  36°,  and  ray-based 
procedures  like  WKBJ  and  Cagniard-de  Hoope  are  tedious  to  implement.  Of  course,  the 
primary  advantage  of  building  up  the  large  ray  files  is  that  the  ray-based  methods  are  more 
readily  modified  to  account  for  lateral  heterogeneity  (e.g.,  Helmberger  et  al.  1985). 
Because  of  the  difficulty  of  identifying  and  including  all  significant  rays  in  the  WKBJ 
synthetics,  and  because  our  primary  interest  here  is  to  understand  basic  characteristics  of 
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the  PP  phase,  we  use  the  reflectivity  method  to  compute  the  remaining  synthetic 
seismograms  in  this  study. 

The  crust  exerts  a  strong  influence  on  PP  waveforms  both  by  reflecting  energy 
from  its  underside  and  by  generating  peg-leg  multiples.  Figure  9  illustrates  the  effects  of 
variations  in  crustal  thickness  on  P  and  PP  waveforms.  The  synthetics  are  computed  using 
the  constant  lid  velocity  model  of  Figure  4,  with  a  single  crustal  layer  having  a  thickness  of 
either  20  or  40  km.  Both  long-period  (top  2  rows)  and  broadband  (bottom  2  rows) 
synthetic  seismograms  are  shown  to  emphasize  the  waveform  differences  resulting  from 
variations  in  crustal  thickness.  The  waveforms  and  relative  timing  of  P  and  PP  phases 
change  significantly  as  the  thickne.ss  of  the  cru-stal  layer  decreases  from  40  to  20  km.  'Hie 
PP/P  amplitude  ratio  is  strongly  influenced  by  absolute  velocities  of  the  shallow  crust, 
which  will  of  course  vary  substantailly.  The  thinner  crust  results  in  reduced  PP-P  travel 
times,  a  very  important  feature  which  could  easily  be  misattributed  to  variation  in  average 
lid  velocity.  The  timing  of  the  PP  arrival  can  also  be  confused  by  the  arrival  of  the 
underside  Moho  reflections  that  precede  and  interfere  with  the  PP  phase,  the  nature  of 
which  is  senstive  to  mantle  velocity  gradients  rather  than  average  velocity  values.  In 
addition,  actual  differences  in  the  average  upper  mantle  P  wave  velocity  structure  will  also 
result  in  PP  amplitude  and  travel  time  variations. 

The  emphasis  here  has  been  on  the  early  crustal  reverberations  that  overlap  the  PP 
arrival  and  produce  its  early  coda.  This  is  generally  called  PP  coupled  PL,  and  is  fairly 
coherent  energy  given  the  degree  to  which  layered  models  can  match  details  of  the 
waveform.  The  later  PP  coda  is  known  to  become  increasingly  incoherent,  with  strong 
scattering  contributions  (Neele  &  Snieder  1991).  It  is  difficult  to  assess  how  much  of  an 
isolated  waveform  should  be  interpreted  in  terms  of  any  average  crustal  structure,  and  one 
must  be  rather  conservative.  In  the  real  earth,  lateral  variations  in  crustal  structure  between 
the  source,  PP  bouncepoint,  and  receiver  regions  is  responsible  for  a  breakdown  in  the 
redundancy  of  crustal  reverberations,  which  is  manifested  in  a  ringing  coda.  Hybrid 
reflectivity  procedures  or  laterally  varying  WKBJ  methods  (e.g.,  Helmberger  et  al.  1985) 
can  be  used  to  try  to  quantify  some  of  the.se  effects,  but  it  is  likely  that  a  scattering 
formalism  (e.g.,  Neele  &  Snieder  1991)  is  the  best  choice  when  the  coda  is  found  to  be 
very  complex.  Broadband  array  measurements  are  the  most  promising  approach  to 
quantifying  the  PP  coda  behavior. 

The  many  factors  affecting  PP  travel  times  and  waveforms  make  it  difficult  to 
attribute  amplitude  and  travel  time  effects  to  one  particular  depth  range  in  the  Earth  without 
some  a  priori  information.  In  many  regions  a  priori  information  on  crustal  thickness  is 
known  from  refraction  profiling.  In  this  ca.se  PP-P  airival  times  and  PP  waveforms  can  be 
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modeled  using  complete  synthetics  to  ascenain  the  average  lid  velocity  structure.  If 
broadband  data  are  available,  constraints  on  lid  gradients  may  also  be  obtained.  Next  we 
illustrate  the  effectiveness  of  using  PP  phases  to  determine  upper  mantle  structure  by 
modeling  long-period  P  and  PP  waveforms  traversing  North  America.  Crustal  thicknesses 
beneath  North  American  are  fairly  well  known  and  we  constrain  this  portion  of  our  models 
using  values  compiled  by  Mooney  &  Braile  (1989). 

3  PP  OBSERVATIONS  AND  WAVEFORM  MODELING 

Many  waveform  modeling  studies  of  the  seismic  velocity  structure  in  the  ppper 
mantle  assume  that  raypaths  travel  primarily  through  laterally  homogeneous  material.  For 
long-period  waves  traversing  isolated  geographic  regions,  this  assumption  is  probably 
reasonable  and  the  one-dimensional  Earth  models  obtained  have  provided  satisfactory  fits 
to  observed  waveforms  (e.g,  Burdick  1981;  Grand  &  Helmberger  1984a;  Walck  1984; 
Lyon-Caen  1986).  We  rely  on  this  approximation  to  model  P  and  PP  waveforms 
propagating  through  similar  tectonic  provinces  in  order  to  determine  lateral  variations  in  the 
upper  mantle  P  wave  velocity  structure  beneath  North  America.  The  advantages  of  using 
multiple  bounce  PP  waves  in  conjunction  with  direct  P  observations  to  determine  upper 
mantle  structure  are  the  increase  in  distance  range  sampling  the  upper  mantle  and  reduction 
in  source  area  anomalies  caused  by  local  structure  or  event  mislocation  through  the  use  of 
differential  arrival  times.  Extension  of  the  distance  range  sampling  the  upper  mantle  is 
essential  for  modeling  upper  mantle  structure  beneath  shield  regions  that  are  often  devoid 
of  earthquake  sources.  The  disadvantages  of  modeling  PP  data  are  the  added  complication 
of  this  phase  as  discussed  in  the  previous  section  and  the  relatively  long  propagation  path 
of  PP,  which  increases  the  likelihood  of  traversing  laterally  heterogeneous  material.  We 
try  to  minimize  the  effects  of  lateral  heterogeneity  on  PP  propagation  by  choosing  source- 
receiver  paths  that  remain  primarily  within  one  tectonic  province. 

Our  data  consist  of  long-period  WWSSN  and  CSN  seismograms  for  paths 
traversing  North  America.  Although  we  have  shown  in  the  previous  section  that  long- 
period  waveforms  have  little  sensitivity  to  velocity  gradients  in  the  upper  mantle,  PP-P 
travel  times  and  gross  features  of  the  PP  waveforms  do  contain  information  about  the 
average  velocity  structure  of  the  upper  mantle,  especially  given  the  independent  constraints 
on  the  structure  provided  by  direct  P  modeling.  Figure  10  shows  our  source-receiver 
geometry  and  Table  1  lists  the  source  parameters  of  the  earthquakes  analyzed.  Locations, 
depths  and  focal  mechanisms  for  all  of  the  sources  were  taken  from  the  literature  references 


noted  in  Table  1.  Earthquake  source  time  functions  were  approximated  as  trapezoids  with 
rise  times  and  source  durations  determined  from  the  inversion  of  teleseismic  P  waves. 

3.1  Upper  Mantle  P  Wave  Velocity  Modeb 

The  earthquakes  we  use  are  all  located  in  the  tectonically  active  regions  of  Alaska, 
Mexico,  and  Central  America  and  are  recorded  at  stations  located  in  the  stable  shield, 
platform,  or  tectonically  active  regions  of  North  America  (Figure  10).  From  S  wave 
studies  of  the  upper  mantle  beneath  North  America  (e.g..  Grand  &  Helmbcrger  1 984a; 
Grand  1987;  Wickens  &  Buchbinder  1980)  and  previous  work  on  P  waves  (e.g.,  Burdick 
&  Helmberger  1977;  Burdick  1981;  Dey-Sarkar  &  Wiggins  1976,  LeFevre  &  Helmberger 
1989;  Wiggins  &  Helmberger  1973;  Walck  1984)  we  expect  upper  mantle  P  wave 
velocities  beneath  the  western  U.S.  to  be  slower  than  beneath  the  Canadian  shield; 
however,  little  is  known  about  the  transition  from  one  province  to  the  other.  To  investigate 
lateral  variations  in  P  wave  velocities  beneath  North  America,  we  construct  a  suite  of  upper 
mantle  P  wave  models  ranging  between  the  shield  model  S25  (LeFevre  &  Helmbcrger 
1989)  and  the  Gulf  of  California  model  GCA  (Walck  1984).  We  then  compare  observed  P 
and  PP  waveforms  to  reflectivity  synthetics  computed  for  each  of  the  models  to  determine 
which  best  fits  the  data.  Since  crustal  thickness  in  most  tectonically  active  regions  is 
thinner  than  in  stable  shield  regions  (Mooney  &  Braille,  1989)  we  assign  crustal 
thicknesses  of  20  km  and  40  km  to  our  tectonic  and  shield  models  respectively;  however, 
we  allow  crustal  thickness  to  vary  if  indicated  by  independent  information  or  if  required  by 
the  data. 

The  suite  of  velocity  models  used  in  our  synthetic  calculations  is  shown  in  Figure 
11.  We  have  slightly  modified  velocities  from  the  published  models  S25  and  GCA  such 
that  all  models  have  the  same  velocities  below  a  depth  of  380  km.  The  velocity  model 
between  380  km  and  700  km  was  modified  from  S25  to  best  fit  data  sampling  all  of  the 
tectonic  regions.  Below  700  km  velocities  are  those  of  PREM.  The  two  end-member 
models  have  about  a  6%  variation  in  P  wave  velocity  in  the  upper  200  km  of  the  mantle, 
resulting  in  PP-P  travel  times  that  can  vary  by  as  much  as  10  s.  Figures  12  and  13 
compare  some  of  the  data  (top  row)  to  synthetic  P  and  PP  waveforms  computed  using  the 
reflectivity  method  for  several  of  the  velocity  models  of  Figure  11.  Examples  of 
waveforms  from  paths  traversing  the  Canadian  shield  are  shown  in  Figure  12  wliile  signals 
traversing  the  North  American  platfomi  and  tectonically  active  regions  are  shown  in  Figure 
13.  The  synthetic  that  best  fits  the  data  is  shown  in  the  second  row  of  each  column,  just 
below  the  data  trace.  The  lower  two  rows  contain  synthetics  for  the  velocity  models  most 
similar  to  the  one  that  best  fits  the  data.  All  traces  in  Figures  12  and  13  are  aligned  on  the  P 


wave  arrival.  Seismograms  from  paths  in  close  proximity  to  one  another  and  traversing  the 
same  geographic  province  are  generally  well  fit  by  the  same  velocity  model,  giving  us 
some  confidence  in  the  stability  of  the  waveform  comparisons. 

For  the  nine  paths  traversing  the  northernmost  Canadian  shield,  the  waveform  data 
are  best  fit  by  the  fastest  model,  the  slightly  modified  version  of  S25  (top  row  of  Figure 
12).  Although  an  average  crustal  thickness  of  40  km  is  reasonable  for  the  Canadian  shield 
(Mooney  &  Braille  1989),  the  model  producing  the  next  fastest  travel  times  is  Shld  (Figure 
1 1)  with  a  20  km  thick  crust.  Figure  12  shows  synthetic  seismograms  for  velocity  model 
Shld,  with  both  a  40  and  20  km  thick  crust.  The  PP-P  differential  travel  time  and  PP/P 
amplitude  ratio  produced  by  the  latter  two  models  are  larger  than  observed  at  the  closer 
stations  (FBC  and  OTT  in  Figure  12),  giving  a  worse  fit  to  the  data  than  model  S25.  By  a 
distance  of  49°,  where  energy  turning  in  he  lower  mantle  constitutes  the  first  arrival, 
model  Shld  with  a  20  km  thick  crust  matches  the  PP-P  travel  time  but  provides  an  inferior 
fit  to  the  observed  waveform  when  compared  with  model  S25.  The  waveform  at  this 
distance  recorded  at  station  STJ  from  event  5  contains  a  cl^ar  later  arrival  from  the  AB 
branch  of  the  upper  mantle  triplication  that  is  best  reproduced  by  the  model  S25  synthetic. 

The  bottom  portion  of  Figure  12  shows  data  and  synthetics  for  three  of  the  six 
paths  traversing  the  southern  portion  of  the  Canadian  shield.  These  data  are  best  fit  by 
synthetics  constructed  with  model  Shld  which  has  slightly  lower  velocities  than  S25  in  the 
upper  200  km  of  the  mantle  (Figure  1 1).  Synthetics  computed  for  the  faster  model  S2.5 
have  PP  waveforms  that  begin  sooner  than  indicated  by  the  data  and  synthetics  produced 
with  the  slower  model  Teca  begin  too  late  and  have  too  small  PP/P  amplitude  ratios.  The 
transition  from  upper  mantle  P  wave  model  S25  to  Shld  appears  to  occur  quite  abruptly 
with  the  seismogram  from  event  5  recorded  at  OTT  best  fit  by  the  faster  model  and  the 
seismogram  from  this  same  event  recorded  at  the  nearby  station  WES  best  fit  by  the  slower 
model.  This  is  compatible  with  the  strong  nonh  to  south  gradients  implied  for  shear 
velocity  structure  by  the  travel  time  station  anomalies  of  Wickens  &  Buchbinder  (1980)  and 
the  tomographic  shear  wave  model  of  Grand  (1987). 

Paths  crossing  the  North  American  platform  also  traverse  crust  with  an  avcragi- 
thickness  near  40  km.  Data  from  three  of  the  five  paths  sampling  the  platfonn  are  show  n 
in  Figure  13  (top  row)  together  with  synthetic  waveforms  indicating  that  model  Tcca  gives 
the  best  match  to  the  ob-served  seismogram.s.  Both  the  travel  times  and  waveforms  of  the 
faster  model  Shld  and  slower  model  Tecb  provide  poorer  fits  to  the  data.  Data,  and 
synthetics  for  three  paths  s'ampling  the  tectonic  province  are  shown  at  the  bottom  of  Figure 
13.  Although  tlie  slowest  model,  GCA  appears  to  be  appropriate  for  many  paths  traversing 
the  tectonically  active  regions  of  the  western  U.S.  and  Canada  (Figure  13),  this  province  is 
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the  most  heterogeneous  and  provides  the  least  satisfactoiy  agreement  between  observations 
and  synthetics.  We  have  found  that  adjacent  raypaths  sometimes  have  waveforms  that  are 
best  fit  by  different  velocity  models.  Crustal  thickness  in  western  North  America  can  vary 
between  20  and  50  km  over  fairly  short  distances  (Mooney  &  Braille  1989),  thus, 
variations  in  crustal  thickness  at  the  source,  receiver,  and  PP  bounce-point  may  account  for 
some  of  the  difficulty  in  matching  observations  with  one-dimensional  models.  Strong 
variations  in  upper  mantle  velocity  .structure,  that  are  not  accounted  for  by  the.se  laterally 
homogeneous  models  must  also  contribute  to  the  wavefomi  mismatch.  Even  given  the 
difficulty  of  finding  an  average  upper  mantle  velocity  model  that  fits  all  of  the  data 
sampling  the  tectonically  active  regions  of  North  America,  it  is  clear  that  some  of  the  data 
do  require  a  velocity  model  as  slow  as  GCA  (Figure  13).  This  is  somewhat  surprising 
considering  that  GCA  was  derived  from  data  sampling  the  actively  spreading  oceanic  ridge 
in  the  Gulf  of  California,  the  tectonic  province  where  we  expect  to  find  the  slowest  upper 
mantle  velocities.  This  is  one  more  line  of  evidence  for  the  extreme  nature  of  the  upper 
mantle  velocity  structure  under  the  basin  and  range  province. 

Waveform  modeling  of  P  and  PP  phases  traversing  North  America  indicates  that  a 
spectrum  of  velocity  structures  exist  beneath  North  America.  The  fastest  velocities  are 
located  beneath  the  nonhern  Canadian  shield,  becoming  progressively  slower  beneath  the 
.southern  shield.  North  American  platform,  and  the  tectonically  active  areas  of  the  western 
U.S..  The  Canadian  shield  and  continental  platform  are  spatially  coherent;  PP  data 
sampling  broad  regions  of  these  provinces  are  well  modeled  by  a  single  one-dimensional 
velocity  structure.  Given  the  complexity  and  structural  sensitivity  of  the  phases 
contributing  to  the  overall  PP  waveform  discussed  in  this  paper  it  is  remarkable  that 
homogeneous  structures  can  account  for  so  much  of  the  wavefield  information.  On  the 
other  hand,  waveforms  sampling  the  tectonically  active  regions  of  the  western  U.S.  are  not 
very  well  matched  using  any  single  one-dimensional  model,  indicating  a  higher  degree  of 
lateral  heterogeneity  at  shorter  horizontal  scale  lengths. 

Upper  mantle  P  wave  velocitiy  models  for  Nonh  America  determined  by  several 
investigators  using  different  data  and  techniques  have  provided  su-ong  evidence  for  the 
existence  of  lateral  variations  in  the  P  and  S  wave  velocity  structure  beneath  the  continent. 
Our  study,  which  utilizes  uniform  data  and  methodology  to  investigate  upper  mantle  P 
wave  velocities  beneath  several  different  tectonic  regions,  supports  the  idea  that  strong 
lateral  velocity  gradients  exist  in  the  upper  mantle  beneath  North  America.  The  average 
upper  mantle  velocities  we  obtain  for  the  different  provinces  of  North  America  agree  well 
with  previous  published  velocity  profiles  determined  for  the  same  regions.  PP  waveform.s 
traversing  the  northern  Canadian  shield  require  P  wave  velocities  in  the  upper  mantle  to  be 


13 


as  fast  as  those  of  model  S25  (LeFevre  &  Helmberger  1989)  while  waveforms  crossing  the 
southern  shield  are  better  fit  by  a  velocity  model  having  slower  upper  mantle  velocities, 
similar  to  the  stable  North  American  continetal  model  S8  (Burdick  1981).  No  single 
velocity  model  can  satisfy  all  of  the  data  traversing  the  tectonically  active  regions  of  North 
America;  however,  upper  mantle  velocities  as  slow  as  GCA  (Walck,  1984)  are  required 
along  some  paths.  We  encountered  little  difficulty  matching  P  and  PP  observations 
sampling  all  of  the  geographic  regions  with  velocity  models  constrained  to  be  identical 
below  a  depth  of  380  km,  and  most  of  our  models  vary  only  in  the  upper  250  km.  More 
extensive  path  coverage  is  required  to  explore  the  abruptness  of  lateral  gradients  in  the 
upper  mantle.  Although  the  long-period  PP  waveforms  are  unable  to  constrain  the  exact 
depth  and  magnitude  of  the  upper  mantle  discontinuities,  it  appears  that  large  variations  in 
these  parameters  are  not  required  by  the  data. 

4  CONCLUSIONS 

A  detailed  study  of  the  PP  phase  at  upper  mantle  distances  has  revealed  that 
underside  reflections  from  the  Moho  and  crustal  reverberations  contribute  significantly  to 
PP  waveform  complexity.  The  failure  of  approximate  synthetic  seismogram  techniques, 
such  as  WKBJ  with  only  primary  rays,  to  include  these  arrivals  results  in  inaccurate  PP 
seismograms  and  has  been  responsible  for  unsatisfactory  agreement  between  PP 
observations  and  synthetics.  More  complete  waveform  synthesis  of  PP  signals  providc.s 
both  improved  modeling  of  the  data  and  potential  for  e,xtracting  additional  constraints  on 
mantle  structure.  In  the  presence  of  a  positive  velocity  gradient  beneath  the  Moho, 
whispering  gallery  phases  will  arrive  prior  to  PP  at  distances  less  than  about  34“,  while 
negative  lid  gradients  lead  to  enhanced  underside  Moho  reflections  at  distances  from  34°  to 
38°.  The  short-period  or  broadband  amplitudes  of  these  underside  Moho  reflections  are 
sensitive  to  the  lid  velocity  gradient,  providing  a  valuable  diagnostic  of  upper  mantle  lid 
structure.  Determination  of  the  velocity  gradient  in  the  lid  and  its  lateral  variations  is  critical 
for  understanding  the  composition  of  the  upper  mantle,  and  lid  gradients  also  cxen  a  strong 
influence  on  the  frequency  characteristics  of  regional  phases.  PP-P  differential  travel  times 
and  gross  PP/P  amplitudes  are  mainly  sensitive  to  a  combination  of  crustal  properties  and 
average  upper  mantle  velocity  structure.  Crustal  reflections  and  conversions  at  the  source, 
receiver  and  PP  mid-path  bounce  point  interfere  to  produce  large  amplitude  PP  coda  that  is 
sensitive  to  crustal  thickness.  Given  constraints  on  the  crustal  thickness,  the  PP-P 
differential  travel  time  can  be  used  to  resolve  average  upper  mantle  velocity  smicture. 


Using  long-period  waves  traversing  various  tectonic  portions  of  North  America,  we 
have  shown  that  if  crustal  thickness  is  known  a  priori  the  average  lid  velocity  can  be 
determined.  We  find  that  a  range  in  average  P  wave  velocity  models  are  appropriate  for  the 
upper  mantle  beneath  North  America.  The  highest  velocities  are  located  beneath  the 
northern  Canadian  shield,  with  velocities  becoming  progressively  slower  beneath  the 
southern  Canadian  shield,  the  stable  continental  platform  and  the  tectonically  jictive 
regions.  Data  from  both  the  shield  and  platform  areas  are  well-modeled  with  one¬ 
dimensional  velocity  structures,  indicating  that  broad  tectonically  stable  regions  of  the 
upper  mantle  are  fairly  homogeneous.  Data  traversing  the  tectonically  active  areas  of  North 
America  can  not  be  as  satisfactorily  modeled  with  a  single  velocity  .structure,  but  the 
average  mantle  velocities  are  certainly  lower  than  beneath  the  craton  and  platfomi.  Active 
tectonic  processes  occurring  in  these  regions  appear  to  have  resulted  in  small  scale  vekx:ity 
heterogeneity. 
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Table  1.  Earthquake  Source  Parameters 


Event 

No. 

Date 

Lat 

CM) 

Long 

CW) 

Depth 

(km) 

ntf, 

Strike 

(*) 

Dip 

(•) 

Rake 

0 

tlt2t3 

(s) 

Reference 

1 

5/17/64 

59.40 

142.70 

10 

5.1 

329 

15 

135 

1 1 1 

PJ 

2 

7/1/67 

54.44 

157.94 

12 

6.1 

94 

72 

5 

122 

HJ 

3 

6/15/68 

5.64 

82.60 

5 

6.0 

3 

82 

180 

.SIS 

P 

4 

8/18/70 

60.70 

145.38 

15 

5.8 

348 

15 

164 

222 

PJ 

5 

7/3/73 

57.99 

138.04 

13 

6.1 

295 

24 

84 

222 

PJ 

6 

4/7/83 

8.0 

82.69 

15 

5.5 

311 

22 

113 

1 1 1 

HVD 

7 

5/9/83 

8.28 

82.93 

7 

5.5 

332 

31 

142 

232 

HVD 

9 

10/13/84 

15.06 

94.24 

16 

6.1 

290 

22 

70 

.5  1 .5 

HVD 

10 

12/20/84 

11.44 

86.16 

15 

5.4 

296 

21 

84 

161 

HVD 

locations  fiom  ISC,  magnitudes  from  ISC  (before  1^83)  or  NEIC  (after  1983),  depths  and  source 
time  histcxies  (tl,t2,t3}  from  tvavefonn  modeling,  focal  mechanisms  fiom  P:  Pennington.  1981; 
HJ:  House  and  Jacob,  1983;  PJ:  Perez  and  Jacob,  1980;  HVD:  Harvard  Centroid  Mcnnent  Tensor. 
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Figure  1.  Long-period  P  and  PP  observed  seismograms  (u^  trace).  WKBJ  synthetics 
(middle  trace),  ai^  reflectivity  synthetics  (bonom  trace)  for  the  Octtrfier  18, 1967  Arctic 
earthquake.  TTie  primary  ray  \VKBJ  synthetic  PP  waveforms  computed  for  model  S25 
(LeFcvre  &  Helmberger  1989)  do  not  strongly  resemble  the  data.  Reflectivity  synthetics 
for  model  S25  differ  markedly  from  the  WKBJ  synthetics  and  more  closely  resemble  the 
data.  WWSSN  and  CSN  station  names  and  distances  firom  the  source  are  listed  above  each 
trace.  AB  and  EF  upper  mantle  triplication  branches  are  marked  where  well  separated. 
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Figure  2.  Ccmparison  between  reflectivity  synthetics  (dashed  line)  and  primaiy  ray 
WKBJ  synthetics  (solid  line)  for  P  and  PP  waveftmns  fw  a  shallow  (ISkm)  vertical  dij^ 
slip  source  at  distances  of  42®  and  48®  for  four  upper  mantle  models  schematically  shown 
here.  Note  that  the  synthetics  agree  when  no  crust  is  included  and  the  degree  of  mismatch 
is  greatest  fw  the  mo^l  with  both  a  crust  and  upper  mantle  discontinuities  (disc). 
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Figure  3.  Synthetic  seismograms  computed  for  a  shallow  (18  km)  vertical  dip-slip 
source  and  the  negative  lid  velocity  gradient  model  of  Figure  4  computed  by:  WKBJ 
including  direct  P,  pP,  sP  and  PP  precursors  reflecting  up  to  diree  times  from  the  bottom  of 
the  Moho^  (PmP);  WKBJ  with  direct  P,  pP,  sP,  PP,  pPP  and  sPP  (PP);  WKBJ  with  the 
combination  of  both  raysets  (PmP+PP);  and  reflectivity  (REF).  WKBJ  synthetics 
computed  with  the  more  complete  rayset  (PmP+PP)  more  closely  resemble  the  initial  part  of 
tlie  reflectivity  synthetics. 
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Figure  5.  Schematic  ray  diagram  of  PP  phase  (PP),  underside  PP  Moho  reflection 
turning  below  the  high  velocity  lid  (PmP),  and  whispering  gallery  phase  (Pn,P(WG)) 
hugging  the  Moho  and  turning  in  the  positive  velocity  gradient  in  the  lid. 
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Figures.  Lx)ng-period  reflectivity  synthetics  computed  for  a  shallow  venical  dip  slip 
source  and  the  vel^ity  models  of  Figure  4  having  negative  gradient,  positive  gradient,  and 
constant  lid  velocities.  Whispering  gallery  arrivals  PniP(WG)  (indicated  by  the  arrows) 

precede  PP  and  are  visible  out  to  distances  of  about  34<>  in  the  synthetics  computed  for  the 
positive  lid  gradient  iiuxlels. 
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Figure  7.  Reflectivity  synthetics  as  in  Figure  6  but  convolved  with  a  broadband 
instrument  rcs^nsc.  WIjispering  gallery  phases  PniP(WG)  (indicated  with  arrows)  for  the 
positive  gradient  velocity  model  are  much  more  pronounced  in  the  broadband 
seismograms. 


rijjMrcS.  Synthetic  seismograms  computed  for  a  shallow  (18  km)  vertical  dip-slip  source 
and  a  slight  positive  vcl(x;ity  gradient  in  the  up}Xir  2(K)  km  computed  using:  WKBJ 
including  only  direct  P  and  PP  rays  and  their  associated  depth  phases  (WKBJl);  WKBJ  for 
the  same  rayset  with  additional  rays  turning  in  the  upper  400  km  of  the  mantle  and 
bouncing  once  from  the  underside  of  the  Moho  (including  the  first  whispering  gallery 
modej  as  well  as  all  rays  having  an  extra  reflection  or  conversion  from  both  the  free  surface 
and  the  Moho  (WKBj2);  complete  reflectivity  (REF).  Inclusion  of  the  more  complete 
rayset  (WKBJ2)  causes  the  WKBJ  synthetics  to  more  closely  resemble  reflectivity 
synthetics  which  include  all  reflections  and  conversions  in  a  layered  Earth. 
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Figure  9.  Long-period  (top  two  traces)  and  broadband  (lower  two  traces)  reflectivity 
synthetics  computed  for  a  shallow  vertical  dip-slip  source  and  the  constant  lid  velocity 
model  of  Figure  4  underlying  a  single  crustal  layer  with  a  thickness  of  either  20  or  40  km. 
Due  to  the  large  contribution  of  crustal  multiples  to  the  PP  phase,  changes  in  crustal 
thickness  have  a  profound  effect  on  PP  wavcfomis.  Waveform  variations  resulting  from 
variations  in  crustal  thickness  are  more  clearly  visible  in  the  broadband  synthetics. 
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Figure  11.  Suite  of  upper  mantle  P  wave  velocity  profiles  used  to  model  data  from  North 
American  paths. 
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Figure  12.  Observed  and  synthedc  P  and  PP  seismograms  for  pauis  traversing  the 
northOT  (top  p^el)  and  southern  (bottom  panel)  portions  of  the  Canadian  shield.  The  data 
are  displayed  in  the  top  row  of  each  column  with  the  remaining  rows  indicating  the 
synthetics  that  best  match  the  data  (.second  row)  and  those  computed  with  the  two  velocity 
mndci.s  most  similar  to  the  best  fitting  model  (lower  two  rows).  Waveform  data  sampling 
the  nonhem  portion  of  the  Canadian  Shield  are  best  fit  by  the  fastest  model.  82.5,  while 
data  sampling  the  southern  shield  rc(]uirc  a  .slower  model,  Shid. 
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Figure  13.  Similar  comparisons  to  those  in  Figure  12  indicating  that  P  and  PP 
wavefoims  sampling  the  continental  platform  (top  panel)  require  a  velocity  model  slower 
than  the  southern  shield  region  (Teca)  and  that  data  from  the  tectonically  active  areas  of 
North  America  (bottom  panel)  sample  upper  mantle  as  slow  as  that  found  beneath  the 
ot^anic  ridge  in  the  Gulf  of  California  (model  Gca). 
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